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The generation of quantum entanglement of macroscopic or mesoscopic bodies in mechanical 
motion is generally bounded by the thermal fluctuation exerted by their environments. Here we 
propose a scheme to establish stationary entanglement between two mechanically oscillating mirrors 
of a cavity. It is revealed that, by applying a broadband squeezed laser acting as a squeezed-vacuum 
reservoir to the cavity, a stable entanglement between the mechanical mirrors can be generated. Us¬ 
ing the adiabatic elimination and master equation methods, we analytically find that the generated 
entanglement is essentially determined by the squeezing of the relative momentum of the mechanical 
mirrors, which is transferred from the squeezed reservoir through the cavity. Numerical verification 
indicates that our scheme is within the present experimental state of the art of optomechanics. 
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I. INTRODUCTION 

Quantum entanglement, as a cornerstone in under¬ 
standing many phenomena in quantum world [1-5], 
serves as a necessary resource in various practical ap¬ 
plications of quantum information processing, such as 
quantum algorithms [6], quantum teleportation [7], and 
quantum crytography [8]. In the past decade, tremen¬ 
dous efforts have been devoted to generate entanglement 
in microscopic systems [9-18]. Recently, the generation 
of entanglement in macroscopic and mesoscopic objects 
and the study of quantum mechanical features in these 
scales have attracted much attentions [19, 20]. 

An optomechanical system supplies an ideal platform 
to explore quantum features of macroscopic or meso¬ 
scopic objects in mechanical motion [21]. Advances in 
this field raise a fundamental question: Whether me¬ 
chanical systems in macroscopic scale exhibit quantum 
behavior? People desire to see under what conditions 
it is feasible to generate nonclassical entangled states in 
macromechanical oscillators. It was found that, similar 
to the microscopic system case, entanglement between re¬ 
mote mirrors can also be generated via optical measure¬ 
ment based on the entanglement swapping idea [22-25] 
and via the coherent interactions between the fields in 
two adjacent cavities [26-28]. However, the emergence of 
quantum effects in macroscopic objects is generally be¬ 
lieved to be bounded by the thermal fluctuation. There¬ 
fore, schemes resorting to an efficient precooling to the 
thermal noise have been proposed to establish a stable en¬ 
tanglement between the mirrors of a cavity [29, 30] and 
between the two dielectric membranes suspended inside 
a cavity [31]. Further studies showed that the reservoir 
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engineering technique supplies a nice idea to entangle the 
mechanical systems without resorting to precooling. For 
example, it is found that a stable entanglement between 
the mechanical oscillators in separated cavities [32] and 
between the cavity mirror and atomic ensemble [33] can 
be generated based on the cascade input-output process. 
Other schemes based on engineering the squeezing char¬ 
acters of the reservoirs have also been proposed to entan¬ 
gle the mechanical mirrors in separated cavities [34, 35], 
in ring-cavity [36], and in double-cavity [37] setups. A 
method using multiple-tone coherent driving to the cav¬ 
ity has been used to entangle two mechanical mirrors for 
a single cavity [38, 39]. 

Inspired by these obvious benefits, i.e. the robustness 
to thermal noise and no precooling, of the reservoir engi¬ 
neering in quantum optomechanical control [32-37], we 
in this work propose a scheme to stably entangle the two 
mechanical mirrors of a single Fabry-Perot cavity via the 
squeezed-reservoir engineering. Going from the master 
equation of the whole system and adiabatically eliminat¬ 
ing the degree of freedom of the cavity field, we derive 
a reduced master equation satisfied by the mechanical 
oscillators. Our analytic study on the mechanical entan¬ 
glement quantitatively characterized by logarithmic neg¬ 
ativity reveals that the entanglement generated comes 
from the squeezing of the relative momentum of the two 
mechanical oscillators. A temperature dependent entan¬ 
glement criterion is obtained, which shows the thermal 
fluctuation tolerance of the generated entanglement. The 
physical condition for achieving the maximal entangle¬ 
ment is explicitly obtained from our analysis. 

Our paper is organized as follows. In Sec. II, we show 
the model and derive the reduced master equation by 
adiabatically eliminating the cavity mode. In Sec. Ill, 
the entanglement generation between the mechanical os¬ 
cillators is explicitly studied. The applicability of the 
adiabatic elimination is also verified. In Sec. IV, a sum- 
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FIG. 1. (Color online) Schematic diagram of a Fabry-Perot 
cavity with two mechanically oscillating mirrors in frequen¬ 
cies wi and u> 2 - A coherent laser with frequency cj l and a 
broadband squeezed laser with squeezing degree r around the 
central frequency £Us are injected into the cavity to implement 
the stable entanglement generation between the mechanical 
oscillators. 

mary is given. 


II. SYSTEM AND ADIABATIC ELIMINATION 

We consider a laser in frequency u> l driven cavity with 
two oscillating mirrors (see Fig. 1) [38, 39]. The ideal 
situation without thermal noise of this system has been 
studied in Refs. [40]. The Hamiltonian of the total sys¬ 
tem is (h = 1) 

H = A c c^c+ [wjOjaj+rij(a^+aj)ctc]+ri(c^+c), (1) 

3 = 1,2 

where c and a 3 , respectively, denote the annihilation op¬ 
erators of the cavity field and the two mechanical os¬ 
cillators formed by the mirrors, rjj denote the coupling 
strengths between the cavity field and the oscillators due 
to the radiation pressure, A c = uj c — wl is the detuning 
of the laser freq uency cu l to the cavity-field frequency 
u> c , and = 2 -^/Pk/wl is the driving amplitude with 
the input laser power P and the cavity damping rate k. 
The mechanical oscillators interact with two independent 
reservoirs at same temperature T. The cavity is further 
driven by a broadband squeezed laser with squeezing de¬ 
gree r around the central frequency ws, which acts as a 
squeezed-vacuum reservoir to the cavity field. Then the 
total system is governed by the Born-Markovian master 
equation [41] 

W(t) = -i[H, W(t )] + C m W(t) + C c W{t), (2) 

where W(t) is density matrix of the total system, 
£m- = Ej=i, 2 7j[(»j + • +njV^ Aj -} with n 3 = 

l/[exp (uij/ksT) — 1] and T>a,p = 2 d ■ p — pd ■ — ■ pd 
represents the dissipators of the two mirrors caused 
by their two independent finite-temperature reservoirs, 
C c . = /t[(7V + l)T> e£ i ■ +fVX>et,a • -(Me* 2Ast £> 6)g • +h.c.)], 
with A s = ws — wl, M = coshr sinhr, and N = sinh 2 r 


represents the dissipator of the cavity field caused by the 
squeezed-vacuum reservoir, y, and k are the damping 
rates of the cavity and the oscillators, respectively. From 
Eq. (2), the steady state value of the cavity field and the 
oscillators can be obtained as (c) ss = Vt,/{in — A c ) = a 
and (a,j ) ss = — r)j\a\ 2 /{ojj — ijj). Then Eq. (1) can be lin¬ 
earized into H = Hq + Hj, where Hq = JA =1 i + 

A c etc and H\ = J2j= 1,2 + dj)(ac i ' + a*c). 

In the large damping limit k 7 ^, we can adiabati- 
cally eliminate the degree of freedom of the cavity field 
and get a reduced master equation satisfied by the two 
mechanical oscillator (see Appendix A) 

P(t) = -*[ P(0] + pit) + W k 

3 = 1,2 j,k= 1,2 

x [(&*+& - C 

+(£/* + C )a]p{t)a k - £,J*a k a]p{i) - p(t)a k a]\ 
+(£/* +£ k )ajp{t)a k - £f*a k ajp(t) - £ k p{t)a k a 3 
+(£]“* + £,t)a)p{t)a\ - £,j*a\a]p(t) - t£p{t)a\a], (3) 

where pit) = Tr c [W(f)], = K+l £ ±UJk) + 
with A = A s = A c and F = N\a\ 2 + Ma 2 e 2,lAt . Keeping 
a traceless structure, Eq. (3) preserves the positivity of 
the reduced density matrix p(t). 

The newly emergent third term in the right-hand side 
of Eq. (3) incorporates all the dynamical effects of the 
cavity field on the two mechanical oscillators. It is inter¬ 
esting to see that the cavity field, as a common contact 
“environment”, can not only induce individual dissipa¬ 
tion (with j = k) to each mirror, but also induce inco¬ 
herent interactions (with j ^ k) between the two mirrors 
by the exchange of virtual phonons. Furthermore, be¬ 
sides the thermal dissipation [the second and third lines 
of Eq. (3)], the squeezing-like dissipation [the fourth and 
fifth lines of Eq. (3)] can also be triggered. It is under¬ 
standable based on the existence of the counter-rotating 
terms in H\. In the special case of r = 0 and single oscil¬ 
lating mirror, after dropping the fast rotating squeezing- 
dissipation terms [42] , Eq. (3) reduces exactly to the sim¬ 
ilar form as the one in Ref. [43] . In this case, the rates of 
cooling and heating denoted by the second and third lines 
of Eq. (3) reduce to (£+* + £t )| r=0 oc [n 2 + (A - uy ) 2 ] -1 
and (£F*+fr )| r=0 oc [k 2 + (A + ay) 2 ] -1 , respectively. It 
means that an efficient cooling is realizable when the cool¬ 
ing rate is larger than the heating one by choosing red- 
detuning driving field (i.e. A > 0 ) [43]. When r ^ 0, the 
incoherent interactions and the squeezing effect induced 
by the cavity field can cause the oscillators in squeezed 
state (see Appendix B), which is crucial for generating 
stable entanglement between the two thermally oscillat¬ 
ing mirrors as shown in the following. 

Considering explicitly the scheme configuration in Fig. 
1 , we have the physical condition that the cavity field 
interacts with the two mirrors in a 7r-phase difference 
77 1 = —772 = rjo■ Further assuming the two mirrors are 
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identical, we get w 3 = w 2 = wo, 71 =72 = 70 , and 
hi = h 2 = no- Then Eq. (3) reduces to 

p{t) = -i[H eS ,p(t)\ 

+ 7o[(h 0 + 1 )®a i ,atP(*) + "o®at ia p(*)] 

z -* J ’ J 3 ’ 3 

3=± 

+ho[(£o"* +£o")®a_,a!^(*) ■*" & * + £o 

+(£o"* + Co )^a_,a_p(i) + (Co * + Co")^al,aL^(*)]’ (4) 

where a± = (fii ± a 2 )/\/2 and FT e // = ._ ± Wj-ajoj + 

[*J7§(C 0 - ~ Co"*)®- + h.c.] with w+ = coo and w_ = 
wo — 277Qlm[^g“ + £/“]. It can be seen that the relative 
motion degree of freedom of the two mirrors is decou¬ 
pled to the center-of-mass one and only the former feels 
the existence of the cavity environment. The presence 
of the squeezing Lindblad terms in the last terms of Eq. 
(4) inspires us that a stable squeezing property could be 
induced to the relative motion degree of freedom by the 
cavity environment, via which a stable entanglement is 
hopefully attainable. 


III. ENTANGLEMENT BETWEEN 
MECHANICAL OSCILLATORS 

A. Analytical criterion for the steady-state 
entanglement 

The entanglement of our system can be measured by 
the logarithmic negativity [44] , which is quantified on the 
covariance matrix of X = (xi,pi, x 2 ,p 2 ), 

Vij = (AXiAXj + AXjAXi)/2 (5) 

with AXj = Xj — ( Xj ), Xj = (dj + dj)/\/2, and pj = 
(i aj — dj)/i/2i. The commutation relations [Xi,Xj] = 

and J = ^ q ^ defines the 

symplectic structure of the system, which is further char¬ 
acterized by the symplectic eigenvalues v = (zq,zz 2 ) of 
the matrix *U • V. The Heisenberg’s uncertainty princi¬ 
ple exerts a constraint on zq such that zy ^ 1/2. Thus 
the Peres-Horodecki criterion [45, 46] is rephrased as that 
the state is separable whenever the uncertainty principle 
is still obeyed by the covariance matrix under the partial 
transposition [47], which connects to V as V = A • V • A 
with A = diag( 1,1,1, —1). If the state is entangled, then 
V would violate the uncertainty principle and zy would be 
smaller than 1/2. The logarithmic negativity just quan¬ 
tifies this violation as [44] 

E n = max{0, -log 2 [2min(z?i,z> 2 )]}. (6) 

In general, the analytical relations among the ele¬ 
ments of the obtained V(f) is hard to evaluate. Under 
the condition that two mechanical mirrors are identical, 


iUij with U = 


J 0 
0 J 


we can prove analytically V(t) = ^ ^ j , where 

U 2 (i) = V 21 (t) and C = -VuW ) with 

Uii(t) + Ui 3 (t) = V 22 {t)+V 24 ,{t) = no + 1/2 (see Appendix 
C). Here we have assumed that the oscillators are initially 
in thermal equilibrium with their respective reservoirs. 
Thus there are only three independent variables in V(f). 

Defining V^ 3 ^(t) = ( Vn(t), V 22 (t), V 12 (t) ) T , we can 
calculate from Eq. (4) 


V (3) (t) = M (3) -V (3) (t) + B (3) (t) (7) 


with 


/ 

O 

CM 

1 

m (3) = 

0 2 (2(L - 

\ 

2C- — u> 0 Wo 

B (3) (t) = 

( 0 .4 >+ 


2wo 

7o) 2(2(!_ - w 0 ) 
2(C - 7o) 



,( 8 ) 

(9) 


where + iQ_ 


^^0 l a | 2 _ \ a \ 2 pr _i_ -pi 

k— i(A+a;o) K-\-i(A— ujq) 5 ^ ^ 


770(^1 +£i"*), and <j> = 7 o( 2 no + 1). The initial condition 
is V ( 3 ) (0) = ( n 0 + 1/2, no+ 1/2, 0 ) T . 


Except for solving Eq. (7) numerically to evaluate the 
entanglement, we can also obtain an analytical form of 
Ejy. For this purpose, we need to convert V(f) into a 
normal form V(f) = U(f)-V(f)-U T (f) by a local uni¬ 
tary transformation U(t) = diag(e lB ' Jy / 2 , e lB ° y t 2 ) with 
a v the Pauli matrix. Such transformation leaves the en¬ 
tanglement unchanged. The achieved V(f) corresponds 
to the covariance matrix defined in the rotated quadra¬ 


ture vector X = (xi,pi,X 2 ,p 2 ) with Xj = [dje lB ! 2 + 
h.c.]/\/2 and pj = [dje~ lB ! 2 — h.c.]/\/2 i. Choosing 9 = 


arg[(aiai)(t)], we get a normal form V(t) 


A C 
C T A 


where 


A = diag (Vji(f), V 22 {t)) ,C = diag (Vi 3 (t),U 24 (i))(10) 

under Vu(f) + V\ 3 (f) = U 22 (f) + U 24 (t) = h 0 + 1/2. Then 
the two symplectic eigenvalues of V(f) are 


£>1,2 = {[Vn (t) ± V 13 (t)}[V 22 (t) T U 24 (f)]}U (11) 


In terms of the centgr-of-mass and relative mption 
quadrature operators Q±_ = {x\ ± X2)/V 7 2 and P± = 
(pi ±p 2 )/V 2, we have V n (t) ± V\ 3 (t) = 5Q±{t ) and 
V 22 it) ± V 24 (t) = 6 P±(t). Thus Eq. (11) changes into 

v 1 = [(®o + 1/2 )6Pl(t)]i, (12) 

z>2 = [(h 0 + l/2)5Q 2 _(t)]X (13) 

The identities 5P 2 (t) = 2 V 2 2 (i) — SP 2 and 5Q 2 _(t) = 
2Un {t)-8Q\ lead to SP 2 (t)—5Q 2 _(t) =2 [U 22 (t)-Un(t)]. 
From the definition of V(f), we have V\\(t) — I/ 22 (f) = 
2Re[(aiai)(t)e _i6) ]. Remembering 9 = arg[(5iai)(t)], we 




4 



FIG. 2. (Color online) Entanglement dynamics in different 
squeezing parameter r when A/ojo = 1 (a) and in different 
A when r = 1 (b), where the black dot-dashed lines show 
the results from the steady-state solution. The steady-state 
entanglement En (blue solid line) and the squared variance 
SPl (red dashed line) in different squeezing parameter r when 
T = 0 K (c) and in different temperature T when r = 1 (d), 
where the frequency detuning is chosen as A = cuo and the 
green dotted lines show the value l/[2(2ho + 1)]. P = 4/rW 
and other parameters are given in the main text. 

readily obtain Vii(t) — V 22 (t) = 2 |(aiai)(f)| > 0 . There¬ 
fore, we always have £7 < r> 2 - Thus the non-zero entan¬ 
glement in Eq. ( 6 ) is achievable only when 

i>i < 1/2 =>■ 5Pl(t) < l/[2(2n 0 + 1)]. (14) 

Equation (14) indicates that the entanglement between 
the oscillators is only determined by the variances of their 

relative momentum operator P_(f), which is only sensi¬ 
tive to the environmental temperature. 

From the commutation relation [Q_,P_] = *, we have 
their uncertainty relation SQlSPl > 1/4. Therefore, the 
squeezing property is present when either SQ ‘1 < 1/2 or 
SPl < 1/2. Thus we can conclude from the analytical 
result in Eq. (14) that it is the squeezing in the relative 

momentum P_ that is responsible for the entanglement 
generation of the mechanical oscillators. At zero temper¬ 
ature, i.e. ho = 0 , the entanglement between the two os¬ 
cillators can be established whenever the squeezing shows 
up in P_; while at finite temperature, the generation of 
entanglement needs stronger squeezing in P_. 

B. Numerical verification 

To verify our analytical expectation, we numerically 
calculate the covariance matrix from Eq. (7) and evaluate 
its entanglement. The parameters are chosen as follows: 
The cavity-field frequency is u c = 2i r x 6.98 x 10° Hz, 


its damping rate is k = 27r x 6.2 x 10 6 Hz, the frequency 
of the mechanical oscillators is wo = 27 t x 32.1 x 10 6 Hz, 
their damping rate is 70 = 15 x 10~ 5 k, and the cou¬ 
pling strength is 770 = 27r x 39 Hz. These parameters are 
accessible with the recent experiment [48]. Figure 2(a) 
shows the entanglement evolution in different squeezing 
parameter r. Obviously, when r = 0, which corresponds 
to the broadband driving field being a vacuum reservoir, 
no entanglement can be stimulated. However, with the 
switching on the squeezing of the broadband driving field, 
a notable entanglement can be induced asymptotically 
in a wide range of r. This demonstrates well that the 
squeezing in the broadband driving field to the cavity 
field is responsible for the entanglement generation of 
the two mechanical oscillators. Figure 2(b) shows the 
entanglement dynamics in different detuning A of the 
broadband driving field to the coherent driving field. We 
can see that the best performance of the squeezing field 
on generating entanglement is achieved when the cen¬ 
tral frequency detuning A of the squeezed field matches 
with the frequency of the oscillators. We also can ob¬ 
serve that the entanglement does not tend to a constant 
value in the long-time limit, but a lossless oscillation with 
tiny amplitude. This oscillation in frequency A origi¬ 
nates from the time-dependent dissipator introduced by 
the broadband squeezed field, as shown by £ c in Eq. (Al) 
and in Eq. (7) with the explicit form of Bl 3 )(t) in Eq. 
(Cll). In the following, we take t = Zn/2 A, where Z 
is an arbitrarily large integer to ensure the system arriv¬ 
ing at its steady state, to study the steady-state prop¬ 
erties of the system. Figure 2(c) plots the calculated 
steady-state quantum entanglement En and the squared 
variance SPl{ 00) = 2V^ 2 (°o) — (n 0 + 1 / 2 ) at zero tem¬ 
perature (i.e. no = 0 ) in different squeezing parameter 
r. It indicates clearly that the region where the sta¬ 
ble entanglement is formed matches well with the region 
where SPl < 1/2. At finite temperature, Fig. 2(d) shows 
that, whenever SPl < l/[2(2fio + 1)] is satisfied, a non¬ 
zero entanglement can be established. All these results 
prove the validity of the entanglement criterion (14). It is 
also impressive to find from Fig. 2(d) that the dramatic 
entanglement can be generated even when the environ¬ 
mental temperature is in the order of magnitude of mK, 
where one generally believes that the dominate thermal 
fluctuation would destroy quantum effects. Here it is re¬ 
markable that the quantum correlation is triggered alive 
even in this temperature region. 

Both of our analytical and numerical results reveal that 
the cavity field acts as a quantum bus to transfer the 
squeezing character from the broadband driving field to 
the relative motion degree of freedom of the two mechan¬ 
ical oscillators such that a stable nonclassical correlation 
can be established between them. It suggests an inter¬ 
esting way to generate stable entanglement between the 
mechanical mirrors by engineering the squeezing prop¬ 
erty of the reservoir felt by the cavity field. The phys¬ 
ical reason why the relative-motion operator instead of 
center-of-mass one is squeezed by the cavity field can be 
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FIG. 3. (Color online) (a) The squared variance SF >2 _ as a 
function of r when P = 0.01 (red dotted line), 0.1 (green 
dashed line), and 2/rW (blue solid line), (b) The optimal 
squeezing r opt (black solid line) to obtain the smallest SP 2 (oo) 
and the corresponding steady-state entanglement (pur¬ 

ple dashed line) as a function of the driving power P. The 
black circles are calculated from numerical solution and the 
black solid line is calculated from Eq. (16). Parameters as 
T = OK, A/cuo = 1 have been used and the others are shown 
in the main text. 


understood in the following way. In our system configu¬ 
ration, 771 = —772 = 770 , under which Eq. (1) is converted 
to 

H = A c c^c + wo ci^aj + V2rj 0 Q_c l c + fl(c^ + c). (15) 

i=± 

It indicates that only the relative coordinate of the os¬ 
cillators feels the existence of the cavity field, while the 
center-of-mass coordinate as a dark mode is immune to 
the cavity field. Thus the squeezing character of the 
reservoir is transferred to relative quadrature operators 
via its interaction with the cavity field (see Fig. 6). This 
again justifies the validity of our analytic result in (14). 
The mediation role of the common cavity field here re¬ 
lates to the scheme in Ref. [39] , where an effective squeez¬ 
ing of the two mechanical modes is achieved by applying 
two coherent driving fields to the cavity, while in ours it 
is transferred from the broadband squeezed reservoir via 
the cavity field. 

Since the entanglement originates from the squeezing 
in P_ , a largest entanglement is achievable when SP 2 ( 00 ) 
has a smallest value. To get the best performance of 
our scheme on entanglement generation, we now explore 
the optimal condition of the system parameters on en¬ 
tanglement generation. Figure 3(a) shows the calculated 
5P_( 00 ) as a function of r in different driving power P, 
which indicates that the optimal value of r to get the 
smallest SP_( 00 ) is dependent on the driving power P. 
Via dS 2 P 2 (oo)/dr = 2dV2i{oo) / dr = 0, we have the op¬ 
timal r to make S 2 P 2 smallest (see Appendix C2) 


r °p t _ 


1 _ ] r © • Re[(2*AI — M^ 3 )) -1 ■ B 2 e 2iAt ] 1 

-arctanhj ; 


(16) 

Figure 3(b) shows r opt obtained from numerical calcula¬ 
tion and from Eq. (16) and the corresponding optimal 
steady-state entanglement as a function of the pumping 
power P. As we can see, with the increase of P, r opt de¬ 
creases and the optimal entanglement generated increases 




FIG. 4. (Color online) SPl(oo) evaluated from Eq. (2) with 
(blue solid line) and from Eq. (4) without the adiabatic elim¬ 
ination (red dashed line) when P = 4^W in (a) and 70 /re = 1 
in (b). Parameters T = OK, A/wo = 1, and r = 1 have been 
used and the others are shown in the main text. 


and saturates at a moderate pumping power. This indi¬ 
cates that a moderate pumping suffices the generation of 
a maximal entanglement between mechanical oscillators. 

Before closing our discussion, we verify the validity of 
the adiabatic elimination used in the derivation of Eq. 
(4). Figure 4 plots SP 2 calculated from Eqs. (2) and 
(4) , which corresponds to the cases without and with the 
adiabatic elimination, respectively. It shows that the adi¬ 
abatic elimination works well in a wide parameter regime, 
e.g., at large cavity decay rate and small and moderate 
pumping power. This result proves the validity of our 
above calculation, where the parameters used are within 
the permissible scope of the adiabatic elimination. 


IV. CONCLUSIONS 

In summary, we have proposed a scheme to generate 
stable entanglement between the two mechanically oscil¬ 
lating mirrors of a cavity by engineering the squeezing 
character of the reservoir felt by the cavity. Via adia- 
batically eliminating the degree of freedom of the cavity 
field, a reduced master equation satisfied by the two me¬ 
chanical oscillators is derived microscopically. From this 
master equation, we have analytically found that the gen¬ 
erated entanglement of the two mirrors originates from 
the squeezing of the relative momentum of the two mir¬ 
rors. Our result reveals that the cavity field acts as a 
quantum bus to transfer the squeezing character of its 
reservoir to the relative momentum of the two mechan¬ 
ical oscillators such that a stable entanglement is estab¬ 
lished in their steady state. The numerical verification 
indicates that our proposal is realizable with the present 
experimental technique of cavity optomechanics. 
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Appendix A: The derivation of the reduced master 
equation 

Working in the interaction picture and further intro¬ 
ducing a squeezing transformation, we can recast Eq. (2) 
into [49, 50] 



W{t) = W{t)} + C m W{t) + C vac W(t), (Al) 

where W(t) = Sexp(iH 0 t)W(t) exp(—iH 0 t)S^ with S = 
exp[r(c 2 - c t 2 )/ 2 ], £ vac - = AZ\ a p, H\{t) = A^{t)c + 
h.c. with A(t) = + h.c.] and = 

rjj(a*y/Ne~ lAt +ay/N + le zAt ). Note that A c = A s = A 
has been used in the above transformation. 

Governed by the dissipator £ v ac, the cavity field 
rapidly approaches the steady state (|0)(0|)_ c in the large 
damping limit (k 3> 7 j). It means that W(t) approxi¬ 
mately factorizes as W(t) ^2 Tr c [W(f)] ® (|0)(0|) c . See¬ 
ing the cavity field as “a reservoir”, we can adiabatically 
eliminate its degree of freedom and obtain a reduced 
master equation satisfied by the two oscillators. Explic¬ 
itly, in the dissipation picture W(t) = exp[—£ vac t]W(f), 

Eq. (Al) can be recast into W(t) = ]A(f) + £ m ]IU(t), 

where C\{t)- = -i exp(-£ vac tJ[Hi(t), ■] exp(£ vac t). Un¬ 
der the Born approximation W(t) = /5(f) ® (|0)(0|) c and 
the Markovian approximation, we obtain 

pOO A A 

pit) = £ m p(t) + Tr c / d,TCi(t)Ci(t-T)p(t)(\0)(0\) c , 
Jo 

^ (A2) 

where p(t) = Tr c [W(t)] and C\(t)- = —i[A+(t)C-(t) + 
A-(t)C+(t ) — h.c.]- with A+(t)- = A^(t)-, A-(t)- = 
A(t)-, C-(t )• = exp(—£ vac t)(c-) exp(£ vac t), and C+(t)- = 
exp(—£ V aci)(c^-) exp(£ vac t). Making a time derivative to 
C_(t), we have dC-(t)/dt = — e - ^ vac *[£ V ac,c-]e^ vac *. One 
can easily check [£ vac ,c-] = nc-. Thus 


FIG. 5. (Color online) Dynamical evolution of the mean 
phonon number calculated with (solid lines) and without 
(dashed lines) the adiabatic elimination in different values of 
r. Parameters T = 2.5mK, A/wo = 1, and 70 /ft = 1.5 x 1(U 3 
have been used and the others are shown in the main text. 


which also results in 


C.\{t) = e Kt (-c) + (e~ Kt - e Kt )(c-). (A7) 


From the obtained forms of Eqs. (A3, A4, A 6 , A7), we 
have the non-zero correlation functions of the cavity field 
as the following 

(C-{t)C + (t - r)) = (C + (t)Cl(t - r)) 

= <Cl(t)C|(t - r)> = (CUt)C + (t - r)> = e~* r (A 8 ) 

where (•) = Tr c [-(10) (01) c ]. 

Substituting Eq. (A 8 ) into Eq. (A2), we can obtain 

poo ^ ^ poo 

Tr c / dr£i(t)£i(f-r)p(t)(|0)(0|) c = / dre~ KT 
Jo Jo 

x [A(t)p(t)A\t — t) — A\t)A(t — r)p(i) + h.c.] (A9) 


Remembering the form of A(t) and returning back to the 
Schrodinger picture, we have 



e - KT e- i6ot A(t)p{t)A\t - r)e iflot dT 


= ^2vjVk[{dj +d])p{t){£+d\+ 4 fife)}, (A10) 

j,k 


C-(t) = e _Kt (c-), (A3) 

which also induces 

Cl(t) = e- ret (-c t ). (A4) 

With the similar manner, we have dC+(t)/dt = 
— e~ Cmct [£ va c, cJ-]e Cva - ct . From the commutation relation 
[£ vac , fiF] =2k ■ c' — it can be recast into 

dC+(t)/dt = — 2 «cl(t) + KC_|_(t). (A5) 

In the form of Eq. (A4), we can obtain 

C+{t) = e K ‘(ct-) + (e~ Kt - e Kt )(-c t ), (A 6 ) 


With £ = «±i(£+A) + and F = + 

Ma 2 e 2lAt . In obtaining Eq. (A10), the integral identity 
f Q exp[— hvT — fwr]dr = 1 /(k + ico) has been used. The 
other terms in Eq. (A9) can be calculated in the similar 
manner. Then the final form the reduced master equation 
(3) can be obtained. 

Appendix B: Dynamical evolution of the system 

The dynamics of the system can be studied readily 
by the derived reduced master equation (4), from which 
we can see clearly that the center-of-mass motion is de¬ 
coupled to the cavity mode, while the relative motion 
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FIG. 6. (Color online) Fluctuation dynamics of SP 2 in (a) 
and 5Q 2 _ in (b) calculated with (blue solid lines) and without 
(red dashed lines) the adiabatic elimination. Insets are results 
in the long-time regime. Parameter r — 1 has been used and 
the others are the same in Fig. 5. 


dissipator introduced by the broadband squeezed field. 
Furthermore, we find that the adiabatic elimination per¬ 
forms well during time evolution except for small devi¬ 
ation in the short-time scale. Figure 6 plots the fluctu¬ 
ation dynamics of the relative motion quadrature oper¬ 
ators P- = (a_ — a*_)/\/2i and Q- = (a_ + at)/\/2- 
The squeezing is achievable when either SP 2 (t) < 1/2 or 
SQ^_(t) < 1/2. The results indicate that the quadrature 
operators can be squeezed in long-time limit, which orig¬ 
inates from the transfer of squeezing properties from the 
reservoir. We also observe that behaviors of the quadra¬ 
ture operators are antiphase with each other and the os¬ 
cillations keeps in a frequency 2A. 


is strongly affected by the optomechanical coupling. For 
the relative motion, not only the thermal dissipation, but 
also the squeezing-like dissipation is triggered. To verify 
the validity of Eq. (4), we plot in Fig. 5 the evolution 
of the mean phonon number n pn = Tr[fit a,jp(t)](j = 1, 2) 
in different values of r obtained from Eq. (2) and from 
Eq. (4). We find that the steady-state phonon num¬ 
ber increases with the increase of r. This is due to the 
strengthen of the radiation pressure acting on the mir¬ 
rors, which comes from the increase of the photon num¬ 
ber in the cavity by the squeezed-vacuum reservoir. In 
addition, faint oscillations appear in the long-time limit 
when r / 0, which originates from the time-dependent 


Appendix C: Covariance matrix of the mechanical 
oscillators 

1. The covariance matrix 

Defining a vector \d F ) (t) = (Vu (i),V 22 (*),V 33 (*),V 44 (*) 
yi2(t),V 1 3(t),V 1 4(t),V23(t),V24(t),V34(t)) T and under the 
conditions wi = w 2 = w 0 , 71 = 72 = 70 , rji = ~V 2 = Vo , 
and fii = U 2 = no, we can derive its time evolution 
equation from Eq. (4) 

V (F) (f) = M • V (F) (t) + B(f ) (Cl) 

with 



-270 

0 

0 

0 

2wo 

0 

0 

0 

0 


0 


0 

2 (Cl - 7o) 

0 

0 

2(C1 — wo) 

0 

0 

-2C1 

—2C1 


0 


0 

0 

-270 

0 

0 

0 

0 

0 

0 


2wo 


0 

0 

0 

2(Cl - 70 ) 

0 

0 

-2C 

0 

-2C1 

2(0 

1 - Wo) 

ci 

- U} 0 

w 0 

0 

0 

1 

1 

to 

0 

-Cl 

-Cl 

0 
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0 

0 

0 
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0 
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-Cl 

0 

0 

0 
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0 
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-Cl 

0 

0 c 
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Cl - 270 
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-Cl 


0 

-Cl 

0 

-Cl 

-Cl 

0 

Cl - w 0 

Pi 

2(Cl - 7o) 


-Cl 


0 

0 0 

1 - W 0 

Wo 

0 

-Cl 

0 

-Cl 
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Cl 

- 270 

it) 

= ( 

C> + 2C r , c/>, 

0 + 2f 

i ci 0, - 

e, -e, -2^, 
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where Cl + 


2r ll \ a \ 2 _ \a \ 2 p r . 

K—i(A-\-LJo) o>o ) 5 ^ ^ 


Vo(£i + £i~*)i and = 7o(2no + 1). Considering that 
the mechanical oscillators are initially in thermal states 


with the same temperature as their reservoirs, we have 
Va( 0) = n 0 + \ and %(0) = 0 for ip j. 

First, according to Eq. (Cl), we have 


J 





( V n (t) - V 33 (t) \ 

( 

-270 

0 

2 wo 

0 \ 

d 

V22 (t) — V44P) I 


0 

2(Cl - 70) 

2 (Cl - w 0 ) 

-20 

dt 

V 12 (t)-V 34 (t) ~ 

Cl-wo 

Wo 

Cl - 270 

-Cl 


V Vi 4 (t) — V 23 (t)) 

V 

-Cl 

0 

-Cl 

Cl - 270 / 


• V (4) (t). 


From the initial condition \d 4 )(0) = (0,0,0, 0) T , its dy¬ 
namical solution can be solved as V( 4 )(f) = (0,0, 0,0) T , 


which indicates 

V 11 (t) = V 33 (f), V 22 (t) = V 44 (t), 
V 12 (t) = V 34 P), V 14 (f) = V 23 (t). 


(C2) 


(C3) 

(C4) 




























Further defining V^ 2 ^(t) = (Vs(f), V ,(t)) T with V${t) = 
V 12 {t) + Vu(t) and V ,(f) = Vu(t ) + V 13 (t) — V 22 (t) — 
V 24 (f), and according to Eq. (Cl), we have V^ 2 ^(t) = 

^ 4 2 J° ^ • V( 2 )(t). Under the initial condition 

VC)( 0 ) = ( 0 , 0 ) T , its dynamical solution can be obtained 
as \++(t) = ( 0 , 0 ) T , which indicates 


Vu(t) = -V 12 (t), (C5) 

V u (t) + V 13 (t) = V 22 (t) + U 24 (t). (C6) 


In last, defining V 7 (t) = V\\(t) + Ui 3 (t), we have 

V 7 (t) = -2 l 0 V 7 {t) + 70(2 n 0 + 1) (C7) 

under the initial condition VV(0) = 77.0 + 7 - Its solution 
reads V 7 (t) = no + 7 , which is time-independent. Thus 

Vu (t) + V 13 (t) = h 0 + -■ (C 8 ) 

With the relations (C3, C4, C5, C 6 , C 8 ), we have the 
final form of the covariance matrix 



2. The solution of the covariance matrix 


The time-dependent inhomogeneous term (9) can be 
separated into 

B< 3 4 5 6 7 8 9 10 11 12 > (t) = B 0 + AHBi + M(B 2 e 2iAt + c.c.), 


A- 

II 

0 

M 

j , > 4 >C- 

0 + C+ - — - 

a 

11 

0 

2c;, c; ) T , 



a 

to 

II 

C+, Y) - 

+ Ki = 

2 ^M 2 , 


270 


(Cll) 


2^n l Q | 2 


«+<(A^, 0 ) ± «+i(A- Mo ) ■ With this se P aration > the dy¬ 
namical solution of Eq. (7) can be constructed as 


V ( 3 ) (f) = -Y • D -1 (I - e Dt )-Y _ 1 -(B 0 + iVBi) 

—M[Y-(D - 2iAI)" 1 -(Ie 2lAt - e Dt ) Y B 2 
+c.c.] + Y-e Dt -Y- 1 -V< 3 >(0), (C12) 

where D = Y _1 • • Y and I is a unit matrix. The 

solution (C12) asymptotically approach [51] 

V^ 3 \oo) = -[M ^]" 1 • (B 0 + JVBi) 

- M[(M (3) - 2iAI ) _1 • B 2 e 2iAt + c.c.], (C13) 


which is a periodically oscillating function in frequency 
2A. From Eq. (C13), we can easily check 


U 22 (oo) = ©-V ( 3 ) (oo) (C14) 

where © = ( sin 2 (9/2) cos 2 (0/2) — sind ). 
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